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Experimental data for the uniform susceptibility, magnetization and specific heat for the material 
Cu2(C5Hi2N2)2CLt (abbreviated CuHpCl) as a function of temperature and external field are compared 
with those of three different dimerized spin models: alternating spin-chains, spin-ladders and the bilayer 
Heisenberg model. It is shown that because this material consists of weakly coupled spin-dimers, much 
of the data is insensitive to how the dimers are coupled together and what the effective dimensionality of 
the system is. When such a system is tuned to the quantum critical point by application of a field, the 
dimensionality shows up in the power-law dependences of thermodynamic quantities on temperature. We 
discuss the temperature window for such a quantum critical behavior in CuHpCl. 

PACS: 75.10.Jm, 75.40.Cx, 75.40.Mg, 75.50.-y, 75.50.Ee 



I. INTRODUCTION 

Quantum disorder and spin-gap phenomena in insu- 
lating and doped magnetic systems have attracted much 
interest recently. A number of novel materials have been 
synthesized, which have a gap in the spin excitation spec- 
trum and consequently exhibit activated thermodynamic 
behavior at low temperatures. Examples of such sys- 
tems include CaV n 2 „+i [1], (V0 2 )P 2 7 [2], SrCu 2 3 
[3], CuGe0 3 [4], Na 2 Ti 2 Sb 2 [5] and a number of or- 
ganic materials such as Cu 2 (CsHi 2 N 2 ) 2 Cl4 abbreviated 
as CuHpCl [6,7]. While the relationship of spin-gap be- 
havior to high temperature superconductivity remains 
one of the most intriguing problems in condensed mat- 
ter physics [8,9], the insulating materials are of interest 
in their own right providing rich interplay of quantum- 
chemistry, strong quantum fluctuations, frustration and 
dimensionality crossovers. 

The material CuHpCl is particularly interesting for 
a variety of reasons. First, from the point of view of 
spin-half models, it is truely in the strong coupling limit, 
where pairs of spins are strongly coupled to each other 
forming a spin-dimer and these pairs are then weakly 
coupled to the rest of the system. This provides a test- 
ing ground for strong coupling theories of the quantum 
disordered phase. Secondly, the exchange energy scale is 
small enough so that by application of an external mag- 
netic field one can drive the system through a phase tran- 
sition from the quantum disordered to a magnetically or- 
dered phase. Third, there is potentially, a hierarchy of 
energy scales so that there is a temperature window of 
low-dimensional quantum critical behavior, which at the 



lowest temperatures will crossover to 3-dimensional be- 
havior. 

We present here results of finite temperature strong 
coupling expansions around dimerized Hamiltonians. We 
have developed a method combining conventional many- 
body perturbation theory interaction representation with 
cluster expansion techniques to carry out these expan- 
sions by fully automated computer programs. Techni- 
cal details of these calculations will be presented else- 
where [10]. At high temperatures, this approach is re- 
lated to conventional high temperature expansions. But, 
unlike those, the strong coupling expansions show excel- 
lent convergence down to very low temperatures for a 
range of parameters, and thus allow to reliably compare 
the experimental data with model Hamiltonians over the 
entire temperature range. We have calculated the uni- 
form susceptibility, magnetization, internal energy and 
specific heat of 3 different classes of model Hamiltoni- 
ans by this method. These include two quasi-one dimen- 
sional models, the alternating spin-chain and the two-leg 
spin-ladder, and one quasi-two dimensional model, the 
spin- bilayer. These models have been the focus of nu- 
merous theoretical investigations in recent years [11-13], 
however, ours is the first comprehensive numerical study 
of the field and temperature dependent thermodynamics 
of these models. 

The material CuHpCl consists of quasi-lD polymeric 
chains [14], where two spin-half Copper atoms are rel- 
atively close to each other and are physically quite far 
from other pairs of atoms. Based on the linear polymeric 
structure, quasi- ID models of weakly coupled spin-dimers 
have been favored for this system. However, the physical 



1 



distance for pairs of atoms between different polymers 
and those within a given polymer are comparable. Fur- 
thermore, one can identify exchange pathways that lead 
to interaction between spins on different chains as well 
as those in a given chain. Thus, in the absence of any 
first principles calculation, the spin Hamiltonian for this 
system and the effective dimensionality is not apriori ob- 
vious. The majority of the experimental data on the ma- 
terial has been interpreted in terms of one-dimensional 
dimcrized spin models [6,7,15], where the exchange cou- 
pling for the dimers is about f3K, whereas the coupling 
between spins on neighboring dimers is about 3K. How- 
ever, there are two pieces of experimental data which 
are in conflict with this picture. First, on application of 
a strong magnetic field when the spin-gap vanishes, a fi- 
nite temperature transition has been observed at temper- 
atures of order IK, whereas no such finite temperature 
phase transition can exist in a strictly one-dimensional 
system. Secondly, neutron diffraction spectra on powder 
samples do not show any evidence of Van-Hove singular- 
ities expected from the quasi-lD models [7]. 

In this study, we find that the uniform thermodynamic 
quantities for these models in zero field are very insensi- 
tive to the way in which the dimers are coupled together. 
In particular, the susceptibility data for this material is 
well described by cither an alternating chain model or 
a spin-ladder model or even a bilayer model. This may 
appear surprising as the latter model is not apparently 
consistent with the structure of the material. The fact 
that one can nevertheless obtain a good quantitative de- 
scription of the succptibility measurements is a strong ar- 
gument that thermodynamic properties are determined 
by local effects thus rendering them almost insensitive to 
the dimensionality This is easily understood in leading 
order of perturbation theory: in a strongly dimerized sys- 
tem the Curie- Weiss constant 9„, , which is the important 
parameter controlling the high temperature susceptibil- 
ity, and the singlet triplet excitation gap A, which is the 
relevant energy scale controlling the low temperature sus- 
ceptibility, are altered from their non-interacting dimcr 
values in the same combination independent of how the 
dimers are coupled together. Still, the extent to which 
one is able to fit the experimental data by adjusting the 
parameters in the three models is surprizing. The theo- 
retical calculations of the specific heat in the three cases 
are also close to each other, although they deviate from 
what is found experimentally. 

It should be noted that the asymptotic low temper- 
ature behavior of thermodynamic quantities do indeed 
depend on the spin-dispersion along different directions. 
The dimensionality enters directly into the power-law 
prefactors multiplying the activated behavior [16]. This 
study shows that such prefactors are essentially impossi- 
ble to see in strongly gapped systems. 

We argue that a robust way to determine the dimen- 
sionality of the spin system is by tuning it to the quan- 



tum critical point by application of a magnetic field. The 
spin-gap in zero field, together with the g-factors, de- 
termine the critical field. As the spin-gap disappears, 
the thermodynamic quantities develop power-law depen- 
dences on temperature [17]. Given that these materi- 
als are ultimately 3-dimensional, an important question 
is: Is there a temperature window where the low dimen- 
sional quantum critical behavior can be observed. We 
find that at temperatures above IK, where our expan- 
sions still show good convergence, the asymptotic low 
temperature power-laws are difficult to determine in an 
unbiassed manner. However, in a coarse sense, quantum 
critical behavior begins to set in at temperature of order 
5K . By studying consistency with expected behavior in 
one and two dimensions, one can determine the dimen- 
sionality of the systems below this temperature. This 
crossover temperature scale for the onset of power-laws 
is consistent with recent NMR measurements [15]. 

The quasi-lD models that we wish to study are given 
by the Hamiltonian: 

H = J± Sa,i ■ Sb,i (1) 

i 

+ Ja ^ &B,i ■ 

i 

Here, Sa,i and Ss,i represent the two spins in unit cell i 
of the chain. The antiferromagnetic intra dimer coupling 
J± > is much stronger than the exchange constants <// 
and J a between different dimers. We will focus on the 
two limiting cases of this general model: the alternating 
spin chain with (J; = , J a ^ 0) and the spin ladder 
defined by (J; ^ , J a = 0). 

In order to compare with a quasi-2D system, we con- 
sider the bilayer Heisenberg model, with Hamiltonian: 

T~l = J± ^ Sa,i • Sba + J2 ^ [Si ; , • Saj + Ss :i • Sbj] , 

i <i,j> 

(2) 

where < i,j > are nearest neighbours on a square lat- 
tice, and A and B represent spins in the two layers. The 
coupling of the system to an external field is described 
by 

H ext = -g^ B Hj2(S z A , i + S z Bii ) (3) 

i 

We perform numerical investigations of these models 
by applying finite temperature strong coupling expan- 
sions [10]. The expansion parameter A is given by ra- 
tio of the inter dimer to the intra dimcr coupling, i. e. 
A = J a /J±, Ji/J±, J2/J1. for the alternating chain, 
ladder and bilayer model respectively. The coefficients of 
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FIG. 1. Uniform susceptibility \ m zero field P er s Pi n vs - 
temperature T. The lines are partial sums of series for the 
three models discussed in the text. The curves are almost 
indistinguishable and provide an excellent fit to the experi- 
mental data of Hammar et al. 

the expansion are polynomials in the variables J±/k B T, 
1/f = l/(l + cxp[-gfi B H/k B T]+c^p[g^ B H/k B T}) 
and 1/Zq = 1/ (1 + /exp[- Jj_/k B T] ) , where Z is the 
partition function of an isolated dimer. 

Scries for the uniform susceptibility \, the Magnetiza- 
tion M, the internal energy E and the specific heat C 
are calculated complete to order A 8 , for arbitrary tem- 
perature and magnetic field. The series coefficients will 
be presented elsewhere and made available on the web. 

To begin our comparison with the experimental data, 
we consider first the uniform susceptibility in zero field. 
We note that there are some deviations between the sus- 
ceptibility measurements of different groups. We will 
consider here the measurements of Hammar et al. [7] and 
use them to fix exchange parameters within the three 
models. In Fig. 1 we compare the experimental data 
with the models with parameters J± = 13K, J\ = 3.5K 
for the ladder, J± = 13.5K, J a = 5.1K for the alter- 
nating chain and J± — 13K, J 2 — 1.9K for the bilayer. 
The fit was obtained with g — 2.04 as given by Hammar 
et al. [7]. We note that the agreement is excellent for 
all three models. Different models essentially reproduce 
the experimental data much better than the difference in 
different experimental results. 

We note that similar ambiguity has also been encoun- 
tered in studying the system (VC>2)P207 [2], where many 
initial measurements were interpreted in terms of a spin- 
ladder model but more recent neutron scattering mea- 
surements suggest that the alternating chain model is 
more appropriate. 

In Fig. 2 the zero-field specific heat measurements 
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FIG. 2. Specific heat C in zero field per spin vs. tempera- 
ture T. The lines are partial sums of series for the three mod- 
els discussed in the text compared with experimental data of 
Hammar et al. 

[7] are compared with the theoretical models. We note 
that while the shape of the experimental spectra is very 
similar and its activated low temperature behavior and 
peak position are well reproduced by the series, the over- 
all quantitative agreement with theory is missing. This 
maybe due to incomplete background subtractions or 
the presence of non-magnetic phases in the experimen- 
tal samples. 

Before we move on to field dependent measurements, 
we turn to the calculation of the critical field. The triplet 
dispersion has been calculated to high orders by zero tem- 
perature dimer expansions [13,18]. For the parameters 
chosen, the gap is determined to be 10.0 K for the ladder 
model, 10.24 K for the alternating spin chain and 9.1 K 
for the bilayer model. Which taking g = 2.04 translates 
into critical fields of 7.30 T, 7.45 T and 6.65 T respec- 
tively. Given that experimental observation of critical 
fields of 7.2(1) Tesla by Hammar et al. [7] and approxi- 
mately 7.7 Tesla by Chaboussant et al. [6], this argues in 
favor of the one-dimcnsiopnal models. 

In Fig. 3 we compare the Knight shift measurements 
of Chaboussant et al. [15] with the M/H ratio calculated 
for the ladder model. Again the two models were found 
to be hardly distinguishable, down to quite low tempera- 
tures. We keep the exchange constants from the previous 
comparisons, but the vertical scale is arbitrary. Given 
that the zero field data of the two groups disagree some- 
what, the agreement, found here, is quite good. The 
convergence of these expansions breaks down once the 
gap closes but only at fairly low temperatures, when the 
correlation length becomes large. 

In Fig. 4 we show the temperature dependence of the 
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FIG. 3. Knight shift vs. log(T). Partial sums of series for 
the ladder model at various magnetic fields compared with 
experimental data of Chaboussant et ai. 



magnetization for the two models, when the system is 
tuned to the quantum critical point. The difference is 
quite apparent. The low temperature power-law behav- 
ior is known [19] to be M oc T 1 / 2 in d=l and MocTin 
d=2. Although the asymptotic low temperature behav- 
ior may set in only at extremely low temperatures the 
difference between the two curves is quite apparent. In 
ID the magnetization approaches zero with an infinite 
slope, whereas in 2D it appears to do so with a finite 
or zero slope. By checking consistency with linear or 
square-root behavior one can distinguish the two cases 
quite clearly at temperatures as high as 4K. The mate- 
rial CuHpCl undergoes a 3-dimensional phase transition 
around IK, hence, the window for observing quantum 
critical behavior maybe very limited. We note that in 
recent measurements of the nuclear relaxation rates for 
this material Chaboussant et al. observe a fairly sharp 
upturn around 5K [15], which they interpret as evidence 
for quantum critical behavior. This crossover tempera- 
ture scale is consistent with what we find here. 

Assuming the spin-ladder model, a mean-field treat- 
ment [20] of the transition temperature, would lead to 
interchain coupling of the same order of magnitude as 
the coupling along the polymeric chains (of order IK). 
This large interchain coupling would explain the absence 
of Van Hove singularities in the powder neutron diffrac- 
tion. On the other hand, if the stronger couplings make it 
a quasi-2D system, with let us say weaker coupling along 
the polymeric chains, there could be a finite temperature 
Kosterlitz-Thouless phase transition, which could turn 
into a 3D phase transition even with very weak 3D cou- 
plings. 

In conclusion, in this paper we have presented finite 
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FIG. 4. Magnetisation m\\ per spin vs. temperature T for 
the critical field H c as discussed in the text. 

temperature strong coupling expansions for the uniform 
susceptibility, magnetization, internal energy and specific 
heat of a number of dimerized spin models at arbitrary 
temperatures and fields. These calculations should be 
quite useful in experimental determination of exchange 
parameters for a class of magnetic materials. We also 
showed that uniform thermodynamic measurements in 
zero field are not sufficient for determining the nature 
of the spin Hamiltonian and the dimensionality of the 
system. We found that the susceptibility data for the 
material CuHpCl can be equally well fit by a number of 
different models. 

We have also shown that field tuning the system to the 
critical point, may provide a clear way to determine the 
effective dimensionality of the spin system and to study 
the asociated quantum critical phenomena. Direct mea- 
surements of the spin dispersion by neutron diffraction 
on single crystals should shed more light on the nature 
of the coupling constants in this material. 

The series can be obtained on the WWW. The 
access is via http://brahms.physik.uni-bonn.de/~ nor- 
bert / series /series. html 
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